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Section 1 


GENERAL DESCRIPTION 


INTRODUCTION 


The Tektronix 4051R01 MATRIX FUNCTIONS are contained in a Read Only Memory device 
designed to be used with the Tektronix 4051 Graphic System. The MATRIX FUNCTIONS ROM 
Pack, shown in Fig. 1-1, is a small plastic case (ROM Pack) containing matrix firmware 
routines for the 4051 microprocessor. The firmware routines are special sets of instructions, 
which have been permanently stored on memory chips inside the ROM Pack. The MATRIX 
FUNCTIONS add to the capabilities of the Graphic System, without altering the operation of 
the system, or occupying any of the available RAM (Random Access Memory) space. 


The MATRIX FUNCTIONS enable the Graphic System to perform certain functions which are 
essential to matrix algebra. Once the ROM Pack is inserted into a slot in the 4051 backpack or 
plugged into a ROM Expander Unit, the Graphic System is able to perform five matrix 
functions, DET, IDN, INV, MPY, and TRN. Briefly, these functions provide the determinant 
(DET), the identity matrix ({ON), the inverse and solutions to systems of linear equations (INV), 
the matrix product (MPY), and the transpose (TRN). 


Using the MATRIX FUNCTIONS is more efficient and convenient than creating BASIC 
programs to perform matrix functions. BASIC programs used to perform the same functions 
"soak up” RAM (Random Access Memory) space. For instance, a BASIC program designed to 
perform the same function as INV takes up at least 60 lines of BASIC program, and ties up 
about 2K bytes of memory space. By contrast, the INV function provided by the MATRIX 
FUNCTIONS ROM Pack does not occupy RAM, and thus leaves more space available for 
storing BASIC programs and data. 


Functions provided by the MATRIX FUNCTIONS ROM Pack are executed much faster than 
BASIC programs used to perform the same function. For example, the INV function inverts a 

10X10 matrix approximately 8.5 times faster than the BASIC program which uses the same 
algorithm. 


in summary, the MATRIX FUNCTIONS ROM Pack is a Read Only Memory device which allows 
the Graphic System to perform five matrix functions which would otherwise require BASIC 
program subroutines. Using the MATRIX FUNCTIONS to perform these functions saves 
programming time and execution time, and conserves memory space. 
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GENERAL DESCRIPTION 
SPECIFICATIONS 


SPECIFICATIONS 
POWER REQUIREMENTS 


The 4051R01 MATRIX FUNCTIONS draws all necessary power from the 4051 power supplies. 
Connections to the power supplies are made through the backpack on the rear panel of the 
4051 main chassis. The MATRIX FUNCTIONS device must be inserted into a slot in the 
backpack, or into a ROM Expander Unit, before power is applied to the system. 


Voltage Supplies Current Limit 


+5 Vde 278 MA 
+12 Vde 25 mA 
—12 Vde 20 mA 


TEMPERATURE 
Non-operating: —40°C to +65°C 
Operating: +10°C to +40°C 
ALTITUDE 
Non-operating: 50,000 feet maximum 
Operating: 15,000 feet maximum 
HUMIDITY 


95% non-condensing (storage) 
80% non-condensing (operating) 


VIBRATION (NON-OPERATING) 


0.015” DA-10-50-10 


SHOCK (NON-OPERATING) 


1/2 Sine 11 ms duration, 30 G’s 
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4051R01 MATRIX FUNCTIONS 


GENERAL DESCRIPTION 
SPECIFICATIONS 


PHYSICAL DIMENSIONS (INCLUDING EDGE BOARD CONNECTOR) 


Length: 4.662” 
Width: 2.620" 
Depth: 0.875” 


WEIGHT 


8 OZ. 


STANDARD ACCESSORIES 


1 — Operators Manual P/N 070-2127-00 


INSTALLATION INSTRUCTIONS 


1. Flip the 4051 power switch to the OFF position. 


NEMO Mel a 


: CAUTION 


Inserting any device into the backpack when the 4051 power is ON may cause 
memory to be erased. Make sure that important information in the RAM is stored on 
magnetic tape before turning the power OFF and proceeding. 


2. With the power removed from the system, insert the MATRIX FUNCTIONS into a slot as 
shown in Fig. 1-2. Press down, and at the same time gently rock the plasic housing from side to 


side until the ROM Pack edgeboard connector is firmly seated in the receptacle connector. 


3. Flip the 4051 power switch to the ON position. After a few minutes of warm-up, the system IS 
ready for use. 


@ 1-3 


GENERAL DESCRIPTION 
SPECIFICATIONS 
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Section 2 


MATRIX FUNCTIONS 


INTRODUCTION 


This section contains a detailed explanation of each of the five functions provided by the 
MATRIX FUNCTIONS Rom Pack. The explanations are intended for someone who is generally 
familiar with the operation of the Tektronix 4051, but who is not necessarily an expert in matrix 
algebra. Background information about matrices is sometimes included, when it is helpful for 
understanding the matrix functions and how they can be used. 


Examples and sample programs are made as simple as possible, and a brief list of pitfalls to 
watch out for is included at the end of each function description. This section begins with a 
"Eunction Summary” chart (Table 2-1), which may be used as a quick reference once you have 
learned how to work with the functions. The chart is also included as Appendix B, so you may 
remove it from the manual and use it as a separate sheet. 


TABLE 2-1 
FUNCTION SUMMARY 


EXAMPLE EXPLANATION ALGORITHM 


Returns the transpose 

of a matrix. Each row 

TRN 2-3 B=TRN(A) (or column) in BilJ}) = AJ.) 
A becomes 

the corresponding 

row column) 

in B=TRN(A). 

If A is an IXJ matrix, 

B=TRN(A) is a JXI 

matrix. 


FUNCTION 








Creates a matrix whose 

diagona! elements ldJ)=tifl=J 
1(1,1), 1(2,2}, 1(3,3),... 

are 1's, and all other l(IlJ) =Oifl¥AdJ 
elements are O's. 


CALL “IDN”, | 
Or 
IDN 2-8 


AS = "IDN" 
CALL A$,| 
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FUNCTION ‘0. EXAMPLE EXPLANATION 


INV 


2-15 


2-26 


TABLE 2-1 (cont) 


C=AMPY8B 


B=INV(A) 


B=INV(A) 
followed by 


DET 
or 
X=DET 








Returns the matrix pro- 


duct of two arrays. If A 
is an IXJ matrix, and 

B is a JXK matrix, the 
product C is an IXK 
matrix. 


Returns the inverse of 
the square part of the 
parameter matrix, and 
solves systems of 
equations represented 
by the matrix. 


Returns the deter- 
minant of the square 
part of the matrix 
most recently supplied 
to the INV function. 





ALGORITHM 
K 
C(I,K) = © A(IL ) * B(L,K) 
L=1 


The Gauss-Jordan tech- 
nique is performed in 
place, with full row 

and column pivoting. 


The value of the deter- 
minant is the product 
of the pivot elements 
computed during the 
INV process. 
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The TRN Function 


The TRN Function 


Syntax Form: 


{ Line number ] array variable = TRN array variable 


Descriptive Form: 


[ Line number ] target variable = TRN parameter variable 





PURPOSE 


The TRN function returns the transpose of a matrix. 


EXPLANATION 
WHAT THE TRN FUNCTION DOES 


When the TRN function is performed ona matrix, the result is anew matrix found by making the 
columns into rows, or the rows into columns. For instance, when B=TRN(A) is performed, 
each row in A becomes the corresponding column in B. This is equivalent to saying that each 
column in A becomes the corresponding row in B. For example: 


13 5 1 8 
A= -B=TRN(A)= | 3 -2 
B25" <G a 


In this example, A is a 2X3 matrix. When the TRN function is performed on A, the result is the 
3X2 matrix B shown above. 


Notice thatthe first row of Ais the same as the first column of 8, and the second row of Ais the 
same as the second column of B. Likewise, the first column of Ais the same as the first row of B, 
the second column of Ais the same as the second row of B, andso on. Matrix B is the transpose 
of matrix A. 


When the TRN function is performed, the number of rows and columns are reversed. In the 


above example, A is a 2X3 matrix, so B=TRN(A) is a 3X2 matrix. In general terms, when Aisa 
matrix having K rows and N columns, B=TRN(A) is a matrix having N rows and K columns. 
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The TRN Function 


2-4 


If Ais a matrix consisting of one row, the TRN function returns a matrix consisting of one 
column, and vice versa. For instance: 


6 
A= [61 8] o=tawiay=| 1 


8 


Inthis example, Ais a 1X3 matrix, that is, a row containing three elements. B is the transpose of 
A, and is therefore a 3X1 matrix, a column containing three elements. 


HOW TO USE THE TRN FUNCTION 
Dimensioning the Parameter Variable 


lf a matrix is to be used as a parameter for the TRN function, it must be dimensioned ina DIM 
statement, using two subscripts to indicate the number of rows and columns in the matrix. For 
example, ifa5X6 matrix Ais to be supplied to the TRN function, A must be dimensioned in the 
following statement: 


DIM A(5,6) 


Using only one subscript to dimension the parameter variable causes an error to occur when 
the TRN function is performed. For instance, when the parameter for the TRN function is a 
matrix consisting of one row of four elements, the matrix should not be dimensioned ina 
statement such as DIM A(4). The correct way to dimension the matrix is shown below: 


DIM A(1,4) 


Dimensioning the Target Variable 


When the TRN function is performed, the result is anew matrix, which must be assigned to an 
array variable, called the target variable because it serves as a target” for the new matrix. The 
target variable must be dimensioned in a DIM statement, using two subscripts to indicate the 
number of rows and columns. The subscripts used to dimension the target variable must be the 
reverse of the subscripts used to dimension the parameter variable. For instance: 


DIM A(5,6),B(6,5) 


This statement dimensions a 5X6 matrix A which may be used as a parameter for the TRN 
function. Matrix B may be used as the target variable, because it is dimensioned to be a 6X5 
matrix. 
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MATRIX FUNCTIONS 
The TRN Function 


A Sample Program 


The following program illustrates how the TRN function may be used to find the transpose of a 
matrix: 


100 DELETE A,B 
110 DIM A(2,3),B(3,2) 
120 READ A 

130 DATA 1,3,5,8,—2,0 
140 B = TRN (A) 

150 PRINT B 


This program computes the transpose of a 2X3 matrix A, and assigns the result to target 
variable B. The first line of the program deletes variables A and B from memory. This is a 
precautionary measure to make sure that the variables are in an undefined state before line 110 
is executed. (Previously defined numeric variables cannot be dimensioned as array variables 
without first deleting them from memory.) 


Line 110 dimensions the variable A to be a 2X3 matrix. Variable B is to be used as the target 
variabie, and is dimensioned in the same statement to be a 3X2 matrix. 


When line 120 is executed, the numbers contained in line 130 are assigned to matrix A in row 
major order. The result is shown below: 


7 1 3 5 
a 4 


In line 140, TRN(A) is performed and the result is assigned to the target variable. Finally, line 
150 causes the result to appear on the display. The result is shown below: 


1 8 
B = TRN (A) = oe 
5 O 
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Naming the Target Variable 


The transpose of a matrix A does not have the same shape as matrix A, unless A is square (has 
the same number of rows as columns). For this reason, caution must be used when giving the 
same name to the target variabie as to the parameter variable. For example: 


100 DELETE A 

410 DIM A(5,2) 

120 READ A 

130 DATA 1,2,6,—4,3,0,8,1,-2 
140 A = TRN (A) 

150 PRINT A 


This program causes a SHAPE ERROR message to appear on the display. The SHAPE ERROR 
occurs when line 140 is executed, because the statement 140 A=TRN(A) uses the same 
variable name for the target variable and the parameter of the TRN function. When the 
statement is executed, an attempt is made to assign the transpose of A, a 2X5 matrix, to array 
variable A, which is dimensioned in line 110 to be a 5X2 matrix. A SHAPE ERROR results. 


When matrix A is square (has the same number of rows as columns), the transpose of A does 
have the same shape as A. In this case, the same variable name can be used for the target 
variable as for the parameter of the TRN function. For example: 


100 DELETE A 

110 DIM A(3,3) 

120 READ A 

130 DATA 1,2,6,—4,3,0,8,1,—-2 
140 A = TRN (A) 

150 PRINT A 
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The IDN Function 


Identity matrices have special properties which are discussed in “Additional Comments: 
Properties of Identity Matrices” on the following pages. 


The IDN Function and Non-square Matrices 


A matrix does not have to have a square shape to be used as the target for the IDN function. The 
IDN function assigns values to the elements of anon-square array I, just as it did to the square 
matrices in the preceding examples: diagonal elements (elements I(1,1),1(2,2),....,.(K,K) ) are 
assigned the value 1, and all other elements are assigned the value 0. 


For example, when | is dimensioned to be a 3X5 array, performing the IDN function results in 
the following matrix: 


Oo 2 Oo 
Ooo © 


10 0 
{= 0710 

00 1 
Note that the diagonal elements 1(1,1), (2,2), and I(3,3), areali 1's, and the restof the elements 
are 0's. Also, the last two columns of the matrix have 0's assigned to every element. 


When the target variable is dimensioned to be anon-square matrix, the IDN function produces 
a matrix having at least one row or column filled with O's. The rows or columns of O’s are the 
ones which lie outside the square part of the matrix. (The square part of a matrix is the largest 
square portion of the matrix which includes the upper-left corner.} For example, let matrix | be 
dimensioned as a 6X3 array: 


110 DIM 1I(6,3) 


When the IDN function is performed, matrix | is assigned the following values: 


The shaded area shown above is the square part of matrix | because it is the largest square that 
can be blocked off, starting in the upper-left corner. Notice that the square part is the 3X3 
identity matrix, and all elements outside the square part are O's. 
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The IDN Function 


Syntax Form: 


| { IDN” 
[ Line number ] CALL« . 
string vanable 


7 array variable 


Descriptive Form: 


[ Line number ] CALL caliname_, target variable 





PURPOSE 


The !DN function creates a matrix whose elements are 1’s along the diagonal, and 0's 
elsewhere. 


EXPLANATION 
WHAT THE IDN FUNCTION DOES 
Creating an Identity Matrix 


The IDN function generates a matrix whose elements are 1's along the diagonal, and 0's 
elsewhere. For example, the IDN function may be used to generate the following 3X3 matrix: 


1 0 0 
|= 01 0 
00 1 


1(1,1), (2,2), and I(3,3) are the diagonal elements of matrix |, and have been assigned the value 
1. All other elements are 0's. 


When the IDN function is used to generate a square matrix as in the above example, the result is 
called an identity matrix. tn this example, matrix | is the 3X3 identity matrix. 


The result of the IDN function is assigned to an array variable, called the target variable. The 
result of the IDN function always has the same dimensions as the target variable. For instance, 
if lis to be the target for the result of the IDN function, and is dimensioned ina DIMstatementto 
be a 5X5 matrix, performing the IDN function results in the 5X5 identity matrix shown below: 


10000 
0100 0 
= 00100 
00010 
0000 1 
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Only four statements are needed to form the identity matrix and display the result. in line 110, 
the variable | is dimensioned to be a square matrix having four rows and four columns. Matrix | 
is to be used as the target for the result of the IDN function, and must be previously 
dimensioned ina DIM statement, or an error occurs. When the IDN function is performed in line 
120, matrix | becomes the 4X4 identity matrix; 1’s are assigned to the diagonal elements of I, 
and 0’s are assigned to all other elements. When line 130 is executed, the result appears on the 
display as follows: 


1000 
= 0100 
0010 
000 1 


Matrix | is now the 4X4 identity matrix. 


The dimensions of matrix | may be changed in order to make an identity matrix of a different 
size. For example, the following program creates the 6X6 identity matrix: 


100 DELETE | 
110 DIM 1{6,6) 
120 CALL “IDN” I 
130 PRINT “l="51 


Line 110 contains the only difference between this program and the program that was used to 
make the 4X4 identity matrix. The DIM statement in line 110 determines which identity matrix is 
generated by the program. Here, | is dimensioned to be a 6X6 array, so the program makes the 
6X6 identity matrix shown below: 


100000 
01000 0 
|= 001000 
00010 0 
0: 063001. 
000001 


Syntax Forms 


In the sample programs discussed above, the CALL statement takes the following form: 


120 CALL "IDN” I 
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When the target variable has a non-square shape, the result of the IDN function resembles an 
identity matrix (the square partis an identity matrix). But because of the extra rows or columns 
of 0’s, the non-square matrix does not have the properties of an identity matrix. For an 
explanation of the special properties of identity matrices, refer to the topic “Additional 
Comments: Properties of Identity Matrices” on the following pages. 


HOW TO USE THE IDN FUNCTION 
Dimensioning the Target Variable 


Before performing the IDN function, a target variable must be dimensioned ina DIM statement, 
using two subscripts to indicate the number of rows and columns in the result matrix. For 
instance, if variable A is to receive the result of the DIM function, and the function is being used 
to generate the 3X3 identity matrix, A must be dimensioned as follows: 


DIM A(3,3) 


Using only one subscript to dimension the target variable causes an error to occur when the 
IDN function is performed. 


Naming the Target Variable 


Any valid numeric variable name may be used as the target for the result of the IDN function. An 
array can serve as the target for the IDN function without having numeric values assigned to 
each element. The only requirement for the target variable is that it be previously dimensioned 
in a DIM statement, using two subscripts to indicate the number of rows and columns in the 
result matrix. 


An array variable name which has numeric values assigned to each element may also be used 
as the target variable. However, the result of the IDN function repiaces all previously assigned 
values. 


Sample Programs 


The following program illustrates how the IDN function is used to create the 4X4 identity 
matrix: 


100 DELETE | 
110 DIM 1(4,4) 
120 CALL “IDN", | 
130 PRINT “l= “jl 
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Invalid Command Using a semicolon (;) between "IDN” and the target variable name, 

Argument causes an INVALID COMMAND ARGUMENT error message to be 
printed on the display. A comma must appear in the CALL 
statement, not a semicolon. 


Forgetting to dimension the target variable also causes an INVALID 
COMMAND ARGUMENT error to occur when the ION function is 
performed. 


Using only one subscript when dimensioning the target variable also 
causes an INVALID COMMAND ARGUMENT error to occur when 
the IDN function is performed. Two subscripts must be used when 


dimensioning a matrix which consists of asingle row orcolumn. For 
example: 


100 DIM A(1,10) 
110 CALL “IDN" || 


Here, matrix | is dimensioned correctly using two subscripts (1 and 
10). 
ADDITIONAL COMMENTS: PROPERTIES OF IDENTITY MATRICES 


Multiplying by an Identity Matrix 


When a square matrix is multiplied by the identity matrix of the same size, the result is the 
original matrix. That is, the following is true for any square matrix A: 


AMPYI=IMPYA=A 


In this expression, MPY is matrix multiplication and | is the identity matrix having the same 
dimensions as A. 


For example, let A be a 2X2 matrix and | be the 2X2 identity matrix as shown below: 


- Eo oe ba 


Multiplying A by | (or | by A) results in the same matrix A: 


Se En Pe 
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The call name IDN in this statement may be replaced by a string variable, provided the string 
variable has been assigned the value "IDN”. For example: 


100 DELETE | 
110 DIM 1(6,6) 
120 A$ = “IDN" 
130 CALL A$,I 
140 PRINT “I=";] 


This program creates the 6X6 identity matrix, just as the last sample program did. This time, 
however, a string variable name appears in the CALL statement. The string variable A$ is 
assigned the value "IDN" in line 120, before appearing as the call name in line 130. The 


statements... 
120 A$ = “IDN” 
130 CALL A$,I 


are equivalent to the statement... 
120 CALL "IDN" | 
NOTE 
The important thing to remember about using a string variable in the CALL 


statement, is that the value "IDN” must be assigned to the string variable, ora CALL 
NAME INVALID error results when the IDN function is performed. 


THINGS TO WATCH OUT FOR WHEN USING THE !DN FUNCTION 


Syntax Error Not inserting a comma between "|DN” and the array variable in the 
CALL statement, causes a SYNTAX ERROR. 


Likewise, leaving out the quotation marks around the callnameIDN 
in the CALL statement, causes a SYNTAX ERROR message to 
appear on the display. 


Call Name Invalid Attempting to execute the statement CALL A$,| results ina CALL 


NAME INVALID error message, unless the value “IDN” has 
previously been assigned to A$. (See “Syntax Forms”.) 
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The MPY Function 


Syntax Form: 


| Line number ] array variable = array variable MPY array variable 


Descriptive Form: 


[ Line number ] target variable = parameter variable MPY parameter variable 





PURPOSE 


The MPY function returns the matrix product of two arrays. 


EXPLANATION 


WHAT THE MPY FUNCTION DOES 


Performing Matrix Multiplication 


When the MPY function is performed, the result is a new matrix. For example: 
1 2 2 8 
A= == 
los} -8=([34] 


= Wee 2 28]  +[— 8 16 
rae F al E a 7 | sa] 


The elements of the new matrix C are found by multiplying each element of the Ith row of the 
first matrix by the corresponding element of the Jth column of the second matrix, then adding 
the values. The sum is the I,Jth element of the product matrix. 


For instance, to do this matrix mulitiplication by hand, start with |I=1 and J=1 and multiply the 
elements of the first row of A by the corresponding elements of the first column of B, then add: 


1#2 +243 = 8 
This is the element found in the first row, first column of C: 


C(1,1) =8 
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weve (3°) El - Ga) 


(Refer to MPY for how to perform these calculations.) 

The above rule only holds when A is square and | is the identity matrix of the same size as A. If} 
is a non-square matrix generated by performing the IDN function on a non-square target 
variable |}, thenitisnottrue that A MPY |] = |[MPYA = A. Instead, the productofAand | (orof! 


and A) may return a matrix which resembles A, buthas chopped off part of A, or added rows or 
columns of O's to A. 


When an Identity Matrix is the Result of Performing Matrix Muitiplication 


One matrix is the inverse of another matrix if their product results in an identity matrix. In other 
words, matrix B is the inverse of matrix A if the following is true: 


AMPY B=BMPYA=| 


In this expression, MPY is matrix multiplication and | is the identity matrix having the same 
dimensions as A and B. A and B must both be square matrices. 


As an example, let A and B be the 2X2 matrices shown below: 


a= taal 82 13-4] 


The products AMPY Band BMPY A both result in the 2X2 identity matrix. This can be verified 
by performing the MPY function on A and B: 
m4, 
0 1 


4 3 =e oS 
coamrroe [42] [3 


-2 37743] _f10 
D=BMPYA= Feat] E a | =[55] 


(Refer to MPY for an explanation of how to compute these products.) Since both AMPY Band 
B MPY A result in the 2X2 identity matrix, A is the inverse of B. 
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oe 3 24 7 4 3 6 
oie ee | = 2 0 J 
4 
2X3 3X3 


As indicated above, matrix A has three columns and matrix B has three rows. Since the number 
of columns tn A is the same as the number of rows in B, the product C = A MPY B can be 
computed. 


Performing E MPY F results in an error for the following matrices: 


10 3 0 
0 1 0 0 1 0 0 
a 2 02 0 aaa 5 2 | 
0 40 1 a 4 
4X4 3X3 


Here, matrix E has four columns and matrix F has three rows. The product E MPY F does not 
exist, and performing G = E MPY F results in an error, because the number of columns in E is 
not the same as the number of rows in F. 


The fact that A may be multiplied by B does not necessarily mean that B may be multiplied by A. 
For instance, in the firstexample C = A MPY Bmay be performed, butD = B MPY Aresultsin 
anerror. The productB MPY Adoes not exist, because the number of columns in B is not equal 
to the number of rows in A: 


1-2 71 
: oO 2 
ut | 


QO 


3X3 2X3 


Since matrix B has three columns and matrix A has two rows, the product B MPY A does not 
exist. 


As afinalexample, if Ais a2X6 matrix, C = A MPY B may be performed for any matrix B which 
has six rows. Thatis, B may be a6X1 matrix, a6X2 matrix, a6X3 matrix, ora 6X4 matrix, andso 


on. 


Insummary, C= A MPY Bcan be performed whenever B has as many rows as A has columns. 
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Next, make |=1 and J=2. Multiply each element in the first row of A by the corresponding 
element in the second column of B, and add: 


1%8 + 2%4 = 16 
This is the value assigned to the first row, second column of C: 
C(1,2) = 16 
Now, let 1=2 and J=1. Using the second row of A and the first column of B, 
C(2,1) = O#2 + (—3)*3 = —9 
And finally, |=2 and J=2. Using the second row of A and the second column of B, 
C(2,2) = 0*8 + (—3) *4 = —12 
These calculations are summarized in the diagram below. At each step, the shaded row of the 


first matrix and the shaded column of the second matrix are the items used to compute the _ 
shaded element of the product matrix: 


STEP 1: (WAZ \\- (Y | C11) =F; 
STEP 2: HAl\Y Ee |° ” , C12) = YY 
STEP 3: Vwwah Ve »\* Voy "| 


STEP 4: yyy) aA |g » _ C(2.2) = Wy 


= 1#2- 2#3=—8 


AN 


1¥#8- 2%4= 16 


WS 
H 
c 
x 
ine) 
+ 
a 
~ 
* 
GW 
| 
| 
(O 
| 


C(2,1) = Wy 


=0*B8 + (—3)*4=-12 


WS 


When the MPY Function Can Be Performed 


The MPY function can only be performed when the number of columns in the first matrix ts 
equal to the number of rows in the second matrix. For instance, AMPY B can be performed for 
the following matrices: 
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Since B has three columns and A has three rows, the product D = B MPY A can also be 
computed. The.result has as many rows as B and as many columns as A, which means D is a 
2X2 matrix: 


BMPY A=D 


2X3 3X2-2X2 


4 44} 


Notice that both A MPY B and B MPY Acan be performed, but the resulting matrices are not 
the same, and are not even of the same size. 


HOW TO USE THE MPY FUNCTION 


Dimensioning the Parameter Variables 


The MPY function returns the matrix product of two arrays. These parameter arrays must 

previously appear in a DIM statement, using two subscripts to dimension each variable. The 
dimensions of the parameters must be such that the arrays can be multiplied. That is, in order 
to perform A MPY B, the second subscript used to dimension A must be the same as the first 
subscript used to dimension B. For example: 


100 DIM A(7,3),B(3,12) 


Dimensioning the Target Variable 


When the MPY function is performed, a new matrix is generated. The result must be assigned to 
a target variable. The target variable must be dimensioned in a DIM statement, using two 
subscripts to indicate the number of rows and columns in the resulting matrix. In order to 
perform C = A MPY B, the first subscript used to dimension C must be the same as the first 
subscript used to dimension A, and the second subscript used to dimension C must be the 
same as the second subscript used to dimension B. For example: 


100 DIM A(7,3),B(3,12),C(7,12) 
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The Dimensions of the Result 


The resultofA MPY Bis anew matrix which has as many rows as A, and as many columns as B. 
For example: 


3 1-2 
A= Rk ; > , B= E I | ; c=amPyB-|° : 4 
2 0 O = 
2X3 3X3 2X3 


Since A has two rows and B has three columns, there are two rows and three columns in the 
result C. 


Here is another example: 


31 
7 _ f2-1 0 
A= | p= [75 Ol 


7 QO 


3 4 7 2-2 
C=AMPYB= E 1 B ! | _ E 6 | : 
| 7 0 rr 147 0 
a 7 
D=BMPYA= —— | +1/<-= f J 
{ Heo — 16 6 


Since matrix A has two columns and matrix B has two rows, C = A MPY Bcan be computed. 
The result has as many rows as A and as many columns as B, which means that C is a 3X3 = 
matrix: 


AMPYB=C 


3X2 2X3-3X3 
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THINGS TO WATCH OUT FOR WHEN USING THE MPY FUNCTION 


Undefined Variable Attempting to perform the MPY function without assigning a 
numeric value to every element of both parameters, causes an 
UNDEFINED VARIABLE error message to be printed onthe display. 


Syntax Error Not assigning the result of the MPY function to a target variable, 
causes a SYNTAX ERROR: for example, the statement A MPY B is 
nota valid statement. Similarly, the statement A MPY B MPY Cis not 
allowed: the result of A MPY B must be assigned to a target variable 
before attempting to muitiply by C. 


Invalid Function Using the same variable name for the target variable as for one of the 

Argument parameters, causes an INVALID FUNCTION ARGUMENT error 
message to appear on the display. (Refer to the preceding topic 
“Naming the Target Variable”.) 


Likewise, forgetting to dimension the target variable results in an 
INVALID FUNCTION ARGUMENT error when the MPY function ts 
performed. 


Shape Error Attempting to perform the MPY function causes a SHAPE ERROR if 
the parameter matrices cannot be multiplied. Two matrices can only 
be multiplied if the number of columns in the first matrix is the same 
as the number of rows in the second matrix. (Refer to “When the 
MPY Function Can Be Performed”.) 


Incorrectly dimensioning the target variable also causes a SHAPE 
ERROR to occur when the MPY function is performed. When 
performing C = A MPY B, the target variable C must be 
dimensioned to have as many rows as A and as many columns as B. 
(Refer to "Dimensioning the Target Variable”) 


Using only one subscript when dimensioning the parameters or the 


target variable, may cause a SHAPE ERROR when the MPY function 
is performed. 
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A Sample Program 


The following program illustrates how the MPY function can be used to find the matrix product 
of two arrays: 


100 DELETE A,B,C 

110 DIM A(3,3),B(3,2),C(3,2) 

120 READ A,B 

130: DATA. 1,2,0,4,14,2;—1,3,0,5;0,-6,.1,2,0 
140 C=AMPYB 

150 PRINT "A="\A 

160 PRINT “B=";B 

170 PRINT "C = AMPY B=";C 


This program performs matrix multiplication on two arrays A and B, and prints the result. The 
parameter matrices A and B are dimensioned in line 110. Since A is dimensioned to be a 3X3 
matrix, and B is dimensioned to be a 3X2 matrix, the result of A MPY B Is dimensioned in the 
same DIM statement to be a 3X2 matrix called C. When line 140 is executed, A MPY B is 
performed and the result is assigned to array variable C. The next three statements cause the 
results to appear on the display, as indicated below: 


1 2 OQ 5 0 =i 32 
- 4 1 2 B= -6 1 , G=AMPY B= os 1 
=i i3. 2 Q =23° 3 


Naming the Target Variable 


The target for the result of the MPY function must not have the same variable name as either of 
the parameters. That is, the following types of statements are not valid: 


A=AMPYB 
B=AMPYB 


Using the same variable name for the target array as for one of the parameters, results in an 
error. 
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Since the results of A MPY B and B MPY A are usually different, care must be taken when 
determining which of the two products A MPY B or B MPY A is needed. Making the wrong 
choice may result in an incorrect answer, or even an error message: in many cases it is possible 
to perform A MPY B, but not B MPY A because of mismatched dimensions. (Refer to the topic 
“When the MPY Function Can Be Performed” for more information.) 


The Zero Matrix 


One unusual property of matrix multiplication concerns the matrix called 0 because it is filled 
entirely with O's. In ordinary arithmetic, X*#Y = 0 means that either X or Y or both areO. But this 
is not the case in matrix algebra: A MPY B may be 0 without either A or B being 0. Here is an 


example: 
Ps = 1 3 _ {00 
a h aI = E 3 Ea ke 


Here, neither Anor Bis0O, yet the product A MPY B is O. Keep in mind that this can happen when 
using the MPY function, and don't be alarmed when it does. 


The Difference Between * and the MPY Function 


There is more than one way to multiply matrices. One way uses the arithmetic operator *to 
form the element-by-element product of matrices A and B. The other way uses the MPY 

function to form the "matrix product” of A and B. The matrix product and the element-by- 
element product are computed in different ways, and yield different answers. For example: 


1 2 _ 28 
aa bh 3 i Fb fi 
7 1 27 f28 2 16 
y= Re F 3 E 4) F wal 


1 2 28) [8 16 
aa apy B= | 0 a b | 7 3 iD 


Notice that A*B and A MPY B give different results. A¥B is called the element-by-element 
product of A and B, because each element in the result is the product of the corresponding 
elements of A and B. That is, C = A*B is calculated as follows: C(1,1) = A(1,1)*B(1,1) = 
1*#2 = 2, C(1,2) = A(1,2)*B(1 2) = 2*8 = 16, and so on. However, D = A MPY B, the matrix 
product of A and B, Is obtained in a different manner. (Refer to the topic “Performing Matrix 
Muitiplication” for more information.) | 
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Size Error Attempting to perform the MPY function causes a SIZE ERROR 
when one or more of the elements in the result is out of range 
(outside the range +1.0E+308). After a SIZE ERROR message 
appears on the display, the result of performing the MPY function 
may be obtained by entering the target variable name and pressing 
RETURN. Eiements in the result which are out of range (too large or 
too small) contain the largest or smallest number possible, 
+8.988465774E+307. 


Ifa SIZE ERROR occurs during execution of a BASIC program, the 
error may be handled without terminating the program by using the 


ON SIZE THEN ... statement. (Refer to the Graphic System 
Reference Manual for an explanation of ON... THEN... statements.) 


ADDITIONAL COMMENTS ABOUT MATRIX MULTIPLICATION 


A MPY B vs B MPYA 


Matrix muitiplication is not commutative, that is, AMPY B does not generally return the same 
result as B MPY A. For instance: 


A= k 3 a E I | 
c= ampve-(° '] [£3] = [2 2) 


7 [43 O 1] | f6 —11 
cai a oo Fe >| i ” 7 b a 


Notice that matrices C = AMPY Band D = B MPY Aare not the same. This is to be expected 
when performing matrix multiplication. 
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The first matrix in the above expression is called the "coefficient matrix” because its elements 
are the coefficients of x and y in the original set of equations. The second matrix is acolumn 
consisting of the unknowns x and y, and the third matrix is a column consisting of the constants 
which appear to the right in the original set of equations. 


If Ais the coefficient matrix, X is the column of unknowns, and B is the column of constants, the 
above expression can be written more compactly as follows: 


AX=B 


AX is the matrix product of A and X. The expression AX=B is equivalent to the original set of 
equations. This fact may be verified by working out the matrix multiplication as follows: 


— £3 2) [x] _ [3x+2y]7_ 
a= (tt) Gd = Lear 


(Refer to the topic “Performing Matrix Multiplication” for how to compute the product by 
hand.) 


5 
But in this example, B is the column, _ therefore the equation AX = B says the following: 


11 
3x + 2y= § 
and 
4x+ y=11 


These are the original two equations. In other words, the matrix product form is an equivalent 
way to state a system of equations. 


NOTE 


Do not enter a matrix of the type B ] from the keyboard. The elements of a matrix 
must always be assigned numeric values, and attempting to assign variable names 
such as x or y results in an error. The discussion above is provided only to help you 
understand the relationship between linear equations and matrix multiplication. 
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Do not confuse the operation A*B with the matrix product provided by the MPY function. 
Matrix multiplication MPY has a special relationship to systems of linear equations, and has 
many other significant applications. You will find when working with matrices that the words 
“product” and “multiplication” often mean the matrix product and matrix multiplication, as 
provided by the MPY function. in such situations, using the operator *instead of the MPY 
function results in incorrect answers. 


Raising Matrices to a Power 


A matrix may be raised to a power in more than one way. One method uses the operator t or * to 
raise each element of an array to a specified power. Another method finds the powers of a 
matrix by using the MPY function. 


For exampie, At2 creates a matrix whose elements are the square of the corresponding 
elements in A. Likewise, A¥A is the element-by-element product of A with itself, and causes 
each element of A to be squared in the same manner as A!2. However, AMPY A squares A by 
computing the matrix product of A with itself, and gives a different result, as seen in the 
following example: 


QO 1 0 1 2 66 
a * A — =A M — 
n= [22]. ammcaze [21]. cameras [2 3] 


It is important not to confuse these methods of raising matrices to a power. Certain situations 
call forthe "power" of a matrix to be computed by repeatedly performing the MPY function. in 
such cases, using the operator | or * by mistake produces incorrect answers. 


How Matrix Multiplication Relates to Systems of Linear Equations 


A set of simultaneous equations can be represented by a matrix product. For example, 
consider the following equations: 


3x +2y= 5 
4x+ y=11 


An equivalent way to state the equations uses a matrix product: 


p Ly] - Fa 
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The shaded portion in each matrix is called the square part because it is the largest square that 
can be blocked off, starting in the upper left-hand corner of the matrix. 


The shaded portion in B=INV{A) is the inverse of the shaded portion in A: 


WY is the inverse of YW) 


Saying that one matrix is the inverse of another means that when the two are multiplied, the 
resultis an identity matrix. (Refer to the IDN Function topic “Additional Comments: Properties 
of Identity Matrices”.) 


The Extra Columns 


In the previous example, both matrices A and B have two columns which lie outside the square 
part. The extra columns are shown below: 


a=[7 377 / B=INV(A)= & a YY) 


The two extra columns in B represent solutions to two different sets of linear equations. The 
first set of equations is 


2x + gy = 1 

x + 2y=3 
and the first extra column in B represents the solution x = —7 and y = 5. The second set of 
equations is 

2x + 3y=4 

x+2y=0 
and the second extra column in B represents the solution x = 8 and y = —4. 


The coefficients of x and y are the same for both sets of linear equations, and can be found in 
the square part of matrix A: 


Vaso] 


The constants which appear to the right of the equal sign in the first set of equations, appearin 
the first column to the right of the square part in A: 


aA 


Lz 0 
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The INV Function 


Syntax Form: 


[ Line number ] array variable = INV array variable 


Descriptive Form: 


f Line number ] target variable = INV parameter variable 





PURPOSE 


The INV function performs matrix inversion and solves systems of linear equations. 


EXPLANATION 


WHAT THE INV FUNCTION DOES 


The INV function returns a new matrix which has the same size and shape as the original array. 
The square part of the new matrix is the inverse of the square part of the original matrix, and any 
columns which lie outside the square part in the new matrix, represent solutions to sets of 
linear equations. A detailed explanation of what this means is given later. Meanwhile, an 
example is given to illustrate what the INV function does: 


~ £2314 7 =f 2-3-7 8 
A=|i 330) = & 2 saad 


Notice that when the 2X4 matrix A is supplied to the INV function, the result B= INV(A) is also 
a 2X4 matrix. 


The Square Part 


The square part of matrix B is the inverse of the square part of matrix A. The square parts of A 
and 8B are the shaded portions shown below: 


pe Vy: 4 - B=INV(A)= YY : 4 
A 
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The First Set of Equations: 
Parameter Matrix 


The Inverse of the Square Part: 
Solutions to Sets of Linear Equations: 
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The Second Set of Equations: 
Parameter Matrix 
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to Equations 
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Fig. 2-1. An example of How to Interpret the Result of the INV Function. 


Equations 
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Likewise, the constants which appear in the second set of equations, are found in the second 
column to the right of the square part in A: 


23 14% 
A = 
b 23 y 
The solution for the first set of equations is found in the first column to the right of the square 
part in B: 


B = INV(A) = & 4 | 


The solution for the second set of equations is found in the second column to the right of the 
square part in B: 


B77 
B = INV(A) = Ee : a) 


In summary, each extra column in B=INV(A) represents the solution to a set of equations 
which uses the elements of the square part in A for the coefficients of x and y. For every extra 


column in B, there is a corresponding extra column in A which contains the constants that 


appear to the right in the set of equations. Fig. 2-1 summarizes the information obtained from 
this example. 


NOTE 


For those who are accustomed to representing systems of linear equations in matrix 
form, it may be more convenient to express the equations and solutions as follows: 


First System of Equations Solution 


Yad) Le] 


Second System of Equations 


attr] -LoJ 


i I-Es] 
d=] 
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HOW TO USE THE INV FUNCTION 
Dimensioning the Parameter Variable 


The parameter of the iNV function must be a matrix having at least as many columns as rows. 
For example, the INV function can be performed on the following matrices: 


|: f f 2 <3 b 2 ss 
O47 ’ 0-1 4 , 0-1 4 6 


If a matrix is to be used as a parameter for the INV function, it must previously be 
dimensioned in a DIM statement, using two subscripts to indicate the number of rows and 
columns. For example, the statement DIM A{3,4) dimensions a matrix which may be used as 
a parameter for the INV function. 


Since the matrix must have at least as many columns as rows, the second subscript used to 
dimension the parameter must be greater than or equal to the first subscript, or an error will 
occur. For example, if matrix A is dimensioned to be a 4X3 matrix in the statement DIM A(4,3) 
attempting to perform INV(A) results in an error. 


Dimensioning the Target Variable 


When the INV function is performed, a new matrix is generated which has the same size and 
shape as the original matrix. The new matrix must be assigned to a target variable. The target 
variable must be dimensioned ina DIM statement, using two subscripts to indicate the number 
of rows and columns. Since the new matrix has the same size and shape as the original matrix, 
the subscripts used to dimension the target variable must be the same as those used to 
dimension the parameter variable. 


Naming the Target Variable 


The target variable may have the same name as the parameter variable. For example, the 
following statement can be used: 


A = INV(A) 
However, when the statement is executed, the original matrix A is replaced by the new matrix 
generated by the INV function. Care should be taken not to aliow the original matrix A to be 


overwritten in this way, if A is to be used in later calculations. In general, it is good practice to 
use different names for the target variable and the parameter variable. 
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5 CAUTION 


it may be tempting in this example to think of matrix A as representing the following system 
of linear equations: 


2x+3y+z=4 
x+2y+z=90 


This is not a correct interpretation. The INV function always solves systems of N 
equations in N unknowns, where N stands for the number of rows in the parameter 
matrix. Since there are two rows in matrix A above, INV(A) provides solutions to 
systems of two equations in two unknowns (x and y). 


The INV Function and Square Matrices 


When the INV function is performed on asquare matrix, the result is the inverse of the original 
matrix. For example: 


, B=INV{(A) = 





The shaded portions shown above are the square parts of matrices A and B. Notice that the 
square part of A is the entire matrix A, and the square part of B is all of matrix B. 


When the INV function is performed, the square part of the result is the inverse of the square 
part of the original matrix. This means that for this example, matrix Bis the inverse of matrix A. 


When the INV function is performed, extra columns in the result represent solutions to sets of 
linear equations. But in this example, there are no extra columns to the right of the square part 
of A or B, so the result does not provide solutions to sets of linear equations. 


To summarize, when the INV function is performed, the square part of the result is the inverse 
of the square part of the original matrix, and any remaining columns represent solutions to sets 
of linear equations. When the matrix supplied to the INV function has a square shape, there are 
no extra columns, and the resulting matrix is the inverse of the original matrix. 
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The fourth column of matrix B represents the solution for the set of linear equations shown 


below: 
x ea 4 =0 
Ox —y =z = 
ex +y aie = 0 


The coefficients of x, y and z in the equations are the elements found in the square part of A, 
and the constants which appear on the right hand side of the equations are found in the last 
column of A. 


The solution to the set of equations is represented by the last column of matrix B, as follows: 


When the INV Function Can Be Performed 


The Shape of the Parameter 


The INV function can only be performed when the parameter matrix has at least as many 
columns as rows. (Refer to the topic ''Dimensioning the Parameter”.) Attempting to perform 
the INV function on a matrix which has more rows than columns, results tn an error. 


The Square Part of the Parameter 


The INV function can only be successfully performed when the square part of the parameter 
matrix has an inverse. Not every square matrix has an inverse. For instance, the matrix C shown 
below does not have an inverse: 


“(et 


When the INV function is Unable to compute an answer because the square part of the 
parameter matrix has no inverse, an INVALID FUNCTION ARGUMENT error message appears 
on the display. For instance, attempting to perform the INV function results in an error for any 
of the matrices shown below: 


Aid aa). ae 3]. Bes : 5] 
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lf a matrix has been overwritten by the result of the INV function, the original matrix may be 
recovered by performing the INV function again; however, truncation may occur during the 
caiculations, causing the answer to be different from the original matrix. How closely the 
answer resembles the original matrix depends upon the “sensitivity” of the parameter matrix. A 
matrix is called “sensitive” if the values assigned to its elements make the matrix more 
susceptible to the truncation errors that normally occur when arithmetic operations are 
performed. 


Matrices which are very sensitive to truncation error during the INV process are called ill- 
conditioned” with respect to inversion. For discussion of truncation and the condition of a 
matrix, refer to the Appendix at the back of this manual. 


A Sample Program 


The following program illustrates how the INV function may be used to invert a matrix and 
solve a set of linear equations: 


100 DEL A,B 
110 DIM A(3,4),B(3,4) 

120 READ A 

130 DATA 1,0,1,0,2,—1,—2,3,—2,1,1,0 
140 B = INV(A) 


150 PRINT "B = INVA =";B 


This program performs the INV function for the following 3X4 matrix: 


1 0 1 0 
A=|2-1-2 3 
2 1 1 «0 


Line 140 computes INV(A) and assigns the result to the target variable B. When line 150 is 
executed, the result appears on the display: 


i; 4 ¢ 3 
B = INV(A) = 2 3 4 9 
0-1-1 -3 


The square part of the resulting matrix B is the inverse of the square part of the original matrix 


A: 
1 1 1 a 0; i 
2 3 4 is the inverse of 2. 1) SZ 


Oa 2 1 1 
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This program performs the INV function on 3X4 matrices. The statement ON SIZE THEN 220 
alerts the system to the possibility of a SiZE ERROR occurring at execution time. If a SIZE 
ERROR occurs while the program is running, control is transferred to line 220. A message then 
appears on the display, explaining that since the square part of the parameter has no inverse, 
the INV function cannot be performed. In this way, a SIZE ERROR has been avoided, and the 
program continues to ask for more data. 


Solving Systems of Equations 


When the square part of the parameter matrix has no inverse and attempting to perform the INV 
function results in an error, the equations represented by the matrix have no solution. 


For example, the square part (the shaded portion) of the matrix shown below has no inverse, 
and attempting to perform the INV function results in an INVALID FUNCTION ARGUMENT 


c= Wee 


Matrix C represents the following set of two equations in two unknowns: 


xX t2y = 1 
3x +6y =9 


Since no solution exists for the above equations, it makes sense that the INV function results in 
an error. 


As another example, suppose you need to solve the following equations: 


ie 7 Sel 
=% Toe =O 
2X hy 22 "1 


The first step is to represent the equations by a matrix: 


O 1 4 1 
D=]-1 0 3 0 
2 Were 7 


The next step is to perform the INV function on matrix D, in order to obtain a new 3X4 matrix 
whose fourth column represents the solution for the set of equations. 


However, the square part of matrix D has no inverse, and the INV function returns an INVALID 
FUNCTION ARGUMENT error. The INV function is unable to solve the equations. 
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For each of these matrices, the square part (the shaded portion) is the 2X2 matrix C from the 
lastexample. Since C does not have an inverse, attempting to perform the INV function on any 
of the above matrices results in an INVALID FUNCTION ARGUMENT error message. 


When the INV Function Cannot Be Performed 


Avoiding an Error Message 


When the INV function is unable to compute a result because the square part of the parameter 
matrix has no inverse, an error occurs. An INVALID FUNCTION ARGUMENT message 
appears whether the error is caused while executing statements directly, or while running a 
BASIC program. However, if the error results while a program is executing, the BASIC 
interpreter behaves as if a SIZE ERROR has occurred, so that the ON SIZE THEN... statement 
may be used to handle the error without terminating program execution. 


When an ON SIZE THEN... statement is executed in a BASIC program, subsequent SIZE 
ERRORs do not halt program execution; instead, control is diverted to the jine number 
specified in the ON SIZE THEN... statement. (Refer to the Graphic System Reference Manual 
for an explanation of ON...THEN... statements.) 


The following program illustrates how the ON SIZE THEN... statement may be used to avoid a 
SIZE ERROR: 


100 ON SIZE THEN 220 

110 DELETE A,B 

120 DIM A(3,4),B(3,4) 

130 PRINT “DO YOU WISH TO PERFORM THE INV FUNCTION?” 

140 PRINT "IF SO, TYPE "1°; IF NOT, TYPE ’0".” 

150 INPUT T 

160 IF T=0 THEN 250 

170 PRINT “ENTER THE ELEMENTS OF A 3X4 MATRIX, IN ROW ORDER:” 
180 INPUT A 

190 B = INV (A) 

200 PRINT " B = INV (A)=":B 

210 GO TO 130 

220 PRINT "THE VALUE OF THE DETERMINANT |S”; DET 

230 PRINT "THE SQUARE PART OF THE MATRIX HAS NO INVERSE.” 
240 GO TO 130 

250 END 
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The DET Function 


Syntax Form: 
[ Line number ] array variable = INV array variable 


| { Line number |] numeric variable = | DET 





NOTE 


As indicated in the syntax form above, the INV function must be performed in order 
to use the DET function. The value of the determinant is computed during the INV 
process. For this reason, it is important to understand the INV function before 
attempting to use DET. 


PURPOSE 


The DET function returns the value of the determinant. 


EXPLANATION 
THE DETERMINANT OF A SQUARE MATRIX 


The determinant of a matrix is a function of the elements in the matrix. The value of the 
determinant is a numeric constant, and can be computed for any square matrix. 


Computing the Determinant of a Square Matrix by Hand 


NOTE 
The following discussion is provided to illustrate what the determinant of a matrix is, 
by showing one method of computing its value by hand. This is notan explanation of 


the algorithm the DET function uses to compute the value of the determinant. 


The value of the determinant is a numeric constant, and can be computed by hand from the 
elements of the array. For example: 


: ee 
A= T= 
ha , DE 2 
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As in the previous example, it makes sense that the INV function cannot compute an answer, 
because no solution exists for the three equations. 


When a matrix is used to represent a system of linear equations as in the preceding two 
examples, the square part of the matrix is called the “matrix of coefficients” because its 
elements are the coefficients of the variables in the equations. When the INV function is unable 
to compute an answer because the matrix of coefficients has no inverse, the system of linear 
equations has no solution (or does not have a unique solution). 


THINGS TO WATCH OUT FOR WHEN USING THE INV FUNCTION 


Undefined Variable Not assigning a numeric value to every element of the parameter 
matrix, causes an UNDEFINED VARIABLE error to occur when the 
INV function is performed. 


Syntax Error Not assigning the result of the INV function to a target variable, 
causes a SYNTAX ERROR: forexample, the statement INV(A) is not 
a valid statement. 


Invaltd Function 
Argument Forgetting to dimension the target variable results in an INVALID 
FUNCTION ARGUMENT error message. 


Likewise, not dimensioning the target variable to have the same size 
and shape as the parameter variable, causes an INVALID 

FUNCTION ARGUMENT error when the INV function is performed. 
The subscripts used to dimension the target variable must be the 
same as the subscripts used to dimension the parameter variable. 


Attempting to perform the INV function when the square part of the 
parameter matrix has no inverse, generally causes an INVALID 
FUNCTION ARGUMENT error message to appear on the display. 
However, if the error is caused during execution of a program, the 
BASIC interpreter behaves as if a SIZE ERROR has occurred. This 
enables the system to respond to the ON SIZE THEN... statement 
and handle the error without terminating program execution. 


Shape Error Attempting to perform the INV function on a matrix which has more 
rows than columns, causes a SHAPE ERROR: the parameter matrix 
must have at least as many columns as rows. 


Attempting to perform the INV function on a matrix which has more 


than 255 rows or more than 255 columns, also causes a SHAPE 
ERROR. 
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Three numbers are obtained in this way, —9, —6, and —12. The value of the determinant of Ais 
the first number, minus the second number, plus the third number: 


DET = -9 +6 -12 = —-15 


The same method can be used to find the determinant when A is a larger matrix. Each element 
along the first row of Ais multiplied by the determinant of the matrix found by covering up the 
row and column containing that element. The value of the determinant of A is equal to the first 
number, minus the second number, plus the third number, minus the fourth number, and so on. 


The computations are tedious for matrices having more than three rows and three columns. 
For example, take A to be the 5X5 matrix shown below: 


rm © 


On — © © 


DET =—43 


> 

\| 
mn - BOS — 
— —- CC) 
oS _' C- O° Ge 
Own wm & 


The first step in computing the determinant is to multiply A(1,1) by the determinant of the 
matrix found by covering up the first row and first column of A: 


Now to compute the value of the determinant of the 4X4 matrix shaded above, begin by 
multiplying the first element in its first row by the determinant of a 3X3 matrix: 





This process continues until all of the required “subdeterminants” have been reduced toa 
numeric value. Notice that each element in the first row of matrix A must be multiplied by a 4X4 
determinant. whose value must be reduced to anumeric constant by performing the necessary 
expansions on 3X3 matrices, and so on. 
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Fora 2X2 matrix as inthis example, the value of the determinant is the product of the diagonal 
elements in the upper-left and lower-right corners, minus the product of the elements on the 
other diagonal: 


DET = 1*4 — 3%2=>—2 


Computing the determinant of a 3X3 matrix by hand is a little harder. 


For example: 


39. b 2 
A= 6 0 34, DET=~—15 
2-1 Q 


One way to begin the calculation is to “expand” along the first row as follows: multiply each 
element along the first row by the determinant of the matrix found by covering up the row and 
column containing that particular element. For instance, multiply the element A(1,1) by the 
determinant of the matrix found by covering up the first row and first column of A: 


sa a 3%[0*0 — (—1)*3]=—9 
2 Yl ~ 


5 aes 2 eee . 
Notice that the determinant of the 2X2 matrix | | is obtained by the method described in 
0 
the first example, that is O#0 —(—1)#*3=3. 


Next, multiply the element A(1,2) by the determinant of the matrix found by covering up the 
first row and second column of A: 


-3G 2 
oY; 


—1 
Zp 


Finally, multiply A(1,3) by the determinant of the matrix found by covering up the first row and 
third column of A: 






; 1% [6*O — 2%3] = —6 





Vi: 


an 


Ui 
yee 3 _ 2e[6"—1 —240] = —12 
Ui, 0 
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Syntax Forms 


When executing statements directly from the keyboard, the value of the determinant may be 
obtained by entering DET and pressing the RETURN key. 


When the DET function appears in a BASIC program, the value of the determinant may be 
obtained froma statement such as 100 PRINT DET. (Astatement such as 100 DET results inan 


error.) 


In either case, the INV function must be performed on the matrix before performing the DET 
function. 


The value of the determinant may be assigned to a numeric variable in a statement such as X = 
DET or 100 X = DET. (Refer to the topic "Assigning the Result to a Target Variable” on the 
foliowing pages.) 


A Sample Program 


The following program illustrates how the DET function can be used to find the determinant of 
the square part.of a matrix: 


100 DEL A,B 
110 DIM A(2,3),B(2,3) 
120 READ A 

130 DATA 1,1,0, 1,2,—-3 
140 B = INV(A) 

150 PRINT "A =""A: 

160 PRINT “DET=";DET 


This program computes the determinant of the square part (the shaded portion) of the 
following matrix: 


=a -s] 


The INV function is performed on matrix A in line 140, and the DET function is used in line 160. 
When line 160 is executed, the value of the determinant appears on the display: 


DET = 1 


The important thing to keep in mind is that the INV function is always performed first, then the 
DET function. 
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It is easy to see the advantage of using the DET function to compute the determinant of sucha 
matrix; the DET function calculates the value of the determinant ina fraction of the time it takes 
to complete the calculation by hand. 


WHAT THE DET FUNCTION DOES 


The DET function returns the determinant of the square part of a matrix which has been used as 
a parameter for the INV function. For example: 


an jy -3|- B = INV(A) = | | ° a) _ DET=1 


As indicated above, after the INV function is performed on matrix A, the DET function is used to 
compute the determinant of the square part (the shaded portion) of A. 


HOW TO USE THE DET FUNCTION 
When to Perform the DET Function 


The DET function is always performed after the matrix has been supplied to the INV function: in 
other words, the INV function is performed first, then the DET function. 


As anexample, the following sequence of statements may be used to find the determinant ofa 
5X5 matrix: 


DEL A,B 

DIM A(5,5),B(5,5) 
INPUT A 

B = INV(A) 

DET 


Notice that the statement B = INV{A) is executed before the statement DET. 
The DET function does not necessarily have to be performed immediately after the INV 
function. No matter how many BASIC statements have been executed since the last use of the 


INV function, the DET function aiways returns the determinant of the square part of the matrix 
most recently supplied to the INV function. 
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SUPPLEMENTARY INFORMATION 


This Appendix provides supplementary information about the suitability of certain matrices for 
inversion and computation of the determinant. The “condition” of a matrix is a ve"y 
complicated topic, and cannot be fully explained in a few pages. The information contained in 
this Appendix is presented as concisely as possible, and is intended only to help you 
understand how the properties of a particular matrix may affect the results of the INV and DET 
functions. 


MATHEMATICAL INVERTIBILITY VERSUS NUMERICAL 
INVERTIBILITY 


In mathematical theory, a matrix either has an inverse or does not have an inverse. If a matrix 
has no inverse, the value of the determinant is zero, and the matrix equation A MPY X = B has 
no solution (or cannot be uniquely solved). But if a matrix has an inverse, the determinant is 
non-zero, and A MPY X = B may be Solved for any B. 


In practice, however, a matrix is called non-invertible when the sequence of arithmetic 
operations used to compute the inverse reaches a point past which the calculations cannot 
proceed further. This definition is not the same as the mathematical definition of invertibility. 
Many matrices are invertible in the mathematical sense, but not in the computational sense. 
For instance, the following matrix is mathematically invertible, but the inverse cannot 
necessarily be computed if ¢ is much smaller than 1: 


Tee Pare 
aa ¥ on 


Although theoretically an inverse exists for this matrix, most machines are unable to compute 
an answer. 


Conversely, a mathematically non-invertible matrix may be “inverted” on most machines 
without causing an error. For example: 


a 1 
ai2 a 
in theory, this matrix has no inverse, yet in practice most algorithms used to compute inverses 


return an answer if 1/a cannot be represented with complete precision (a=3 is an example of 
this). 


A-1 
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Assigning the Result to a Target Variable 


The DET function always returns the determinant of the square part of the matrix most recently 
supplied to the INV function. This means that as soon as the INV function is performed ona 
new matrix, the result of the DET function changes. 


To store the value of the determinant for later use, the result of the DET function may be 
assigned to a numeric variable in a statement such as X = DET or 100 X = DET. The numeric 
variable is called the target variable because it serves as a “target” for the value of the 
determinant. 


Naming the Target Variable 


itis important to keep in mind that the result of the DET function is a numeric constant, and 
should be assigned to anumeric variable (that is, a variable which has not previously appeared 
ina DIM statement). Thus, the target variable should not have the same name as an array. For 
example, if the INV function is performed on a 3X3 matrix named A, the assignment A= DET 
replaces all nine elements of A by the numeric constant resulting from the DET function. 


THINGS TO WATCH OUT FOR WHEN USING THE DET FUNCTION 


Syntax Error Entering a statement such as 100 DET is a SYNTAX ERROR. A 
correct statement is 100 PRINT DET. (When statements are entered 
without a line number, the keyword DET may be used by itself.) 


Undefined Variable Forgetting to perform the INV function before the DET function, 
causes an UNDEFINED VARIABLE error message to appear on the 
display. The INV function must be performed first, then the DET 
function. 


NOTE 


The DET function always returns the determinant of the square part of the matrix 
most recently supplied to the INV function. Forgetting to perform the INV function 
before the DET function may result in an incorrect answer. For example, when 
finding the determinant of a matrix called B, if DET is performed before INV{B), the 
answer which appears on the display is not the determinant of matrix B; it is the 
determinant of whatever matrix was fast supplied to the INV function. 
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THE CONDITION NUMBER 


The condition of a matrix with respect to inversion can be measured, and the measure is called 
the condition number k. For any matrix C, the condition number may be calculated as follows: 


k =| largest element of C * largest element of INV(C) | 


If k is close to 1, matrix C is well-conditioned with respect to inversion. However, if k is much 
larger than 1, C is less well-conditioned; and if k is extremely large, C is called ili-conditioned., 
and the computed inverse of C is likely to be very far from the true inverse of C. 


The number 1/k is a rough measure of the “distance” from matrix C to the nearest 
mathematically non-invertible matrix. For instance, if the condition number k is 10’, C may be 
made non-invertible by increasing or decreasing at least one of its elements by about 10 * 
times the largest element in the matrix. 


How to Use the Condition Number k to Evaluate the Result of the INV Function 


The condition number may be used as acheck of the accuracy of the result of the INV function. 
If the distance 1/k from C to the closest non-invertible matrix is comparable in magnitude to 
possible errors in the data (the values assigned to the matrix), the result of performing the INV 
function may be so inaccurate that INV(C) is meaningless. For example, if the values assigned 
to the elements of matrix C are physical observations which are only accurate to three places, a 
condition number k of 10° or higher indicates that the result of INV(C) may be too inaccurate to 
be meaningful. Likewise, if the distance 1/k from C to the closest non-invertible matrix is 
comparable in magnitude to the precision of the machine, 10 ", the result of INV(C) may be too 
inaccurate to be of any practical value. This means that if the condition number k is 10°’, the 
distance from C to the closest non-invertible matrix is roughly 10 '*, which is the same as the 
amount of error which might normally occur during any arithmetic operation such as addition 
or subtraction. In general, if the condition number k is 10° or larger, the matrix may be so ill- 
conditioned with respect to inversion that the INV function cannot be expected to compute a 
meaningful result. 


Evaluating the Result of the DET Function 


The value of the determinant is computed during the INV process. For this reason, the result of 
the DET function is affected by the condition of the matrix in the same manner as the INV 
function. When the result of the INV function is inaccurate because the parameter matrix is ill- 
conditioned with respect to inversion, the value returned by the DET function is not likely to be 
close to the true value of the determinant. 


A-3 
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TRUNCATION ERROR 


Truncation is the reason for the difference between what should theoretically happen and what 
actually happens when a machine is used to calculate an inverse. Truncation can occur within 
any machine while arithmetic operations are being performed, and while numbers are being 
converted from decimal to binary representations. Whenever a binary number has more bits 
than the machine can represent, the extra bits are chopped off (truncated). Truncation is 
unavoidable, since there are many numbers which cannot be represented with compiete 
precision on any machine. 


When a sequence of arithmetic operations is performed, some truncation usually occurs. The 
total amount of error depends on the numbers involved in the calculations. For matrix 
inversion, the numbers involved are the vaiues assigned to the elements of the matrix, and the 
“condition” of a matrix is an indication of how these values influence the total amount of 
truncation error. 


THE CONDITION OF A MATRIX 


A matrix C is called "ili-conditioned” with respect to inversion when small changes in the 
elements of C cause large changes in the computed inverse of C, or cause C to be non- 
invertible ina computational sense. It may not be worth attempting to compute the inverse ofa 
very ill-conditioned matrix, because small errors in the origina! data (the values assigned to the 
elements of the matrix) can cause the computed answer to be unrecognizably different from 
the true inverse of the matrix. Ili-conditioned matrix also tend to be more susceptible to 
truncation error than other matrices. 


A matrix which is ill-conditioned with respect to inversion is more susceptible to truncation 
error in the sense that as arithmetic operations proceed to compute the inverse, truncation 
error tends to accumulate faster than it would for better-conditioned matrices. Normal 
amounts of truncation error do "pile up” even for well-conditioned matrices, just because so 
many computations are needed to invert a matrix, and every step in the process provides 
another chance for error to occur. (Thus large matrices are more susceptible to truncation 
error than small ones.) Yet when an ill-conditioned matrix is inverted, much more than the = 
normal amount of truncation error may accumulate, causing the result to be too inaccurate to 
be called an inverse at all. 


The most extreme example of this is anon-invertible matrix, that is, one which has no inverse in 
the mathematical sense. This type of matrix is the worst case of ill-conditioning with respect to 
inversion, and the algorithm used to compute the inverse is able to return an answer, the “0 
answer is completely inaccurate. So much truncation error has accumulated during the 

calculations that the machine appears to have found an inverse, when actually no inverse 

exists. “n 
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Since the size of the condition number indicates that R may not be accurately inverted, itis a 
good idea to check the result by performing R MPY S and S MPY R. How accurate the result 
returned by the INV function is, depends on how closely the elements of the products resemble 
those of the 3X3 identity matrix. 


The products R MPY S and S MPY R are shown below: 


1 0 0 a. =0:5. “2 
M=RMPYS=]0 2 0; , N=SMPYR=j]-6 2 —8 
pe a 8. =e. O55 <2 


Neither product resembles the 3X3 identity matrix, confirming that matrix R is too ill- 
conditioned to be accurately inverted. 


In general, after using the INV function to inverta matrix and solve systems of linear equations, 
or after performing the DET function, it is useful to check the condition number k. An extremely 
large condition number indicates that the parameter matrix is not suitable for matrix inversion 
or computation of the determinant, and that the results of the INV and DET functions should 


not be assumed to be accurate. 
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Thus, if the values assigned to the elements of the matrix are based on physical observations, 
4/k should not be smaller than possible errors in the original data. Likewise, a condition 
number k greater than or equal to 10" indicates that the parameter matrix is so sensitive to 
truncation error that the result of the DET function may be very different from the true value of 
the determinant. 


This means that the value returned by the DET function cannot be used as an indication of 
whether or nota matrix is invertible. In mathematical theory, a matrix has no inverse when the 
value of the determinant is 0; but in practice, if the parameter matrix is non-invertible or ill- 
conditioned with respect to inversion, the value returned by the DET function may not even be 
close to O. 


An Example 


Suppose you want to use the INV function to find the inverse of the following 3x3 matrix: 


1 QO 1 
R= ]|2-1T 5 
3-1 6 


Performing INV(R) results in the following matrix: 


6.251999482E+13 6.254999482E+13 —6.254999482E+13 
S=INV(R) = | —1.876499845E+ 14 —1.876499845E+ 14 1.876499845E+ 14 
—6.254999482E+13 —6§.254999482E+13 6.254999482E+ 13 


Since performing the INV function does not result in an error, you might be tempted to accept 
matrix S as the inverse of matrix R, without checking the accuracy of the result S. 


But a check of the condition number k reveals that matrix R is extremely ill-conditioned with 
respect to inversion: 


| 


=| largest element of R * largest element of S| 
R(3,3)* S(2,3) 
1.125899907 E+15 


li 


l{ 


The condition number k is so large, and R is so ill-conditioned with respect to inversion, that 
the INV function cannot be expected to return an accurate inverse. 
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- FUNCTION SUMMARY 


| eee | EXPLANATION ALGORITHM 


Returns the transpose 
of a matrix. Each row 
(or column) in 

A becomes 

the corresponding 
row column) 

in B=TRN({A). 

If Ais an [XJ matrix, 
B=TRN{A) ts a JXI 
matrix. 


FUNCTION 



















TRN 2-3 B=TRN(A) Bild) = At) 










Creates a matrix whose 
diagonal elements 
1(1,1), 1(2,2), 1(3,3),... 
are 1's, and all other 
elements are O's. 


CALL "IDN", | 
or i= Vitis 


IDN 2-8 
iJ) =Oif lA J 


AS = "IDN" 
CALL A$, 


| 1 | | 





Returns the matrix pro- 
duct of two arrays. If A 
is an IXJ matrix, and 

B is a JXK matrix, the 
product C is an IXK 
matrix. 









K 
C(I,K) = = ACL )*B(L,K) 
L=1 


MPY 2-15 C=AMPYB 


The Gauss-Jordan tech- 














Returns the inverse of 














INV 2-26 B=INV{A) the square part of the | nique is performed in 
parameter matrix, and | place, with full row 
- solves systems of and cojumn pivoting. 
equations represented 
by the matrix. 
B=INV(A) 
followed by Returns the deter- The value of the deter- 
; DET 2-37 minant of the square minant is the product 





of the pivot elements 
computed during the 
INV process. 


DET 
or 
X=DET 


part of the matrix 
most recently supplied 
to the INV function. 
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Appendix C 


ERROR MESSAGES 


CALL NAME INVALID IN IMMEDIATE LINE — MESSAGE NUMBER 32 
— A string variable appears in the CALL statement without being previously assigned 
the value “IDN"”. 
INVALID COMMAND ARGUMENT — MESSAGE NUMBER 12 
— The target variable for the 1DN function has not been dimensioned in a DIM 


statement. 


— The target variable for the IDN function has been dimensioned using only one 
subscript. 


— A-semicolon appears in the CALL statement e.g. CALL “IDN"|| 


INVALID FUNCTION ARGUMENT — MESSAGE NUMBER 27 


— The target variable for the INV, MPY, or TRN function has not been dimensionedina 
DIM statement. 


~~ INVALID FUNCTION ARGUMENT — MESSAGE NUMBER 21 


— The target variable for the INV function is incorrectly dimensioned. 


INVALID FUNCTION ARGUMENT — MESSAGE NUMBER 30 


— The target variable for the MPY function has the same name as one of the 
parameters. 


INVALID FUNCTION ARGUMENT — MESSAGE NUMBER 31 


— The square part of the parameter for the INV function has no inverse. 


SHAPE — ERROR—MESSAGE NUMBER 8 


— The parameter for the INV function has more rows than columns. 


—- The parameter for the INV function has more than 255 rows or more than 259 
columns. 


— The target variable for the MPY or TRN function is incorrectly dimensioned. 


— The target variable or one of the parameters for the MPY function has been 
dimensioned using only one subscript. 
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TEKTRONDG 


4051R05 
BINARY PROGRAM 
LOADER 


Operator’s Manual 


Tektronix, Inc. 

P.O. Box 500 

Beaverton, Oregon 97077 
070-2171-00 


ERROR MESSAGES 


— The parameters for the MPY function cannot be multiplied, because the number of 
columns in the first parameter is not equal to the number of rows itn the second 


parameter. 


— The. parameter for the TRN function has a non-squareshape, and the the target 
variable has the same name as the parameter variable e.g. C=TRN(C). 


SIZE ERROR — MESSAGE NUMBER 1 


— Floating point overflow occurs when the MPY function is performed. 


UNDEFINED VARIABLE — MESSAGE NUMBER 36 


— One or more of the elements of a parameter for the INV, MPY, or TRN function has 
not been assigned a numeric value. 


— The DET function has been performed without first supplying the matrix to the INV 
function. 
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Section 1 


GENERAL DESCRIPTION 


INTRODUCTION 


The Tektronix 4051R05 BINARY PROGRAM LOADER is an accessory to the Tektronix 
4051 Graphic System. Four firmware routines are housed in a Read Only Memory device 
known as a ROM pack. The firmware routines are special sets of instructions for the 

4051 microprocessor. The BINARY PROGRAM LOADER adds to the capabilities of 

the Graphic System without disturbing the normal operation of the system or occupying 


any RAM (Random Access Memory) storage space. 


The BINARY PROGRAM LOADER gives the Graphic System the capabilities to work with 
program files stored in machine dependent binary code. Without the BINARY PROGRAM 
LOADER, the Graphic System can only work with data files stored in machine dependent 
binary code by using the READ and WRITE commands. Now, with the BINARY PRO- 
GRAM LOADER, the Graphic System has four binary program routines—BSAV, BOLD, 
BAPPEN, and LINK. With these routines, you can now save, retrieve, or append a binary 
program file, and you can also use binary program files as sub-routines by using the 


LINK routine. 


When you use the BINARY PROGRAM LOADER, you may not notice anything different 
about your program. Your program is listed on the screen just like it always is. You do 
not see any special code. You may notice that you don’t have to wait as long to store a 
program. The main function of the BINARY PROGRAM LOADER is storage format, 


which is something you never see. 


MACHINE DEPENDENT BINARY CODE 


The term ‘“‘machine dependent binary code” refers to the binary format used by the 
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TEMPERATURE 


Non-operating: 


Operating: 


ALTITUDE 


Non-operating: 


Operating: 


HUMIDITY 


GENERAL DESCRIPTION 


—40°C to +65°C 
+10°C to +40°C 


50,000 feet maximum 


15,000 feet maximum 


95% non-condensing (storage) 


80% non-condensing (operating) 


VIBRATION (NON-OPE RATING) 


0.015" DA-10-50-10 


SHOCK (NON-OPERATING) 


1/2 Sine 11 ms duration, 30 G‘s 


PHYSICAL DIMENSIONS (INCLUDING EDGE BOARD CONNECTOR) 


Length: 4.662” 
Width: 2.620” 
Depth: 0.875" 
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GENERAL DESCRIPTION 


Graphic System to store BASIC programs and data in the internal Random Access Memory 


(RAM}. This binary code is unique to the 4051 operating system. 


Each BASIC command in an ASCII program file must be translated from ASCII code to 
machine dependent binary code in order to be ready for execution. Storing a program 

in machine dependent binary code saves time. For example, the BOLD routine loads 

a program into RAM from a peripheral storage device in approximately 1/4 the time that 
the OLD command takes to load the same program stored in ASCII code. This time is 
saved because the program does not go through a translation. After being loaded, the 


program is ready to be executed or edited just like a program stored in ASCII code. 


From now on, this manual refers to machine dependent binary code as binary code. 


SPECIFICATIONS 


POWER REQUIREMENTS 


The 4051RO05 BINARY PROGRAM LOADER draws all necessary power from the 4051 
power supplies. Connections to the power supplies are made through the backpack on 
the rear panel of the 4051 main chassis or from the ROM Expander Unit. The BINARY 
PROGRAM LOADER ROM pack must be inserted into a slot in the backpack, or into a 
ROM Expander Unit, before power is applied to the system. 


Voltage Supplies Typical Current 
Operating | Non-operating” 

+5 Vde 56 mA 8mA 

+12 Vdc 20 mA OmA 

—12 Vdc 4mA 4mA 


“The ROM pack is inserted but the routines are not being used at this time. 
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Fig. 1-1. 4051R05 BINARY PROGRAM LOADER ROM Pack. 
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2171-1 


GENERAL DESCRIPTION 


WEIGHT 


8 OZ. 


STANDARD ACCESSORIES 


1--Operators’ Manual........ 070-2171-00 


INSTALLATION INSTRUCTIONS (For the 4051 ROM Backpack) 


1. Flip the 4051 power switch to the OFF position. 


CAUTION § 


Inserting any device into the backpack when the 4051 power is ON 





may cause the contents of RAM to be erased. Make sure that impor- 
tant information in RAM is stored on magnetic tape before turning 


the power OFF and proceeding. 


2. With the power removed from the system, insert the BINARY PROGRAM LOADER 
ROM pack into a slot as shown in Fig. 1-2. Press down gently while rocking the 
plastic housing from side to side until the ROM pack edgeboard connector is firmly 


seated in the receptacle connector. 


3. Flip the 4051 power switch to the ON position. After a few moments of warm-up, 


the system is ready to use. 
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BINARY PROGRAM LOADER ROUTINES 


INTRODUCTION 


Each of the four firmware routines provided by the BINARY PROGRAM LOADER is 
explained in this section. If you are already familiar with the operation of the Tektronix 


4051, then these routines are a natural extension of the Graphic System you know. 


By using the Graphic System BASIC language CALL statement, you can specify the 
BINARY PROGRAM LOADER routine you want. The routine name can be specified 
as a string constant in the CALL statement, or the routine name can be assigned to a 
string variable and represented in the CALL statement by that string variable. Refer to 
the 4051 Graphic System Reference Manual for more detailed information on using the 
CALL statement. 


Three of the routines (BSAVE, BOLD, BAPPEN) perform the same basic operation on 
binary program files as their counterpart commands (SAVE, OLD, APPEND) perform 
on ASCII program files. The syntax of the commands and the storage format of the files 
operated on may be different, but the same basic operation is performed. For example, 
the BSAVE routine stores a program on a peripheral device in binary code, while the 
SAVE command stores a program on a peripheral device in ASCII code. The following 


chart shows the relationships. 


ASCIil Files BINARY Files 
SAVE BSAVE 
OLD BOLD 
APPEND BAPPEN 


The explanation of these BINARY PROGRAM LOADER routines includes comparisons 
to the related ASCII file commands. Use the 4051 Graphic System Reference Manual to 


find information on the ASCII file commands. 
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2471-2 


Fig. 1-2. Inserting the BINARY PROGRAM LOADER ROM Pack into a backpack slot. — 
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STORING A PROGRAM ON MAGNETIC TAPE IN BINARY CODE 


To execute the BSAVE routine, the read/write head of the output device must be posi- 
tioned at the beginning of a file marked BINARY or NEW. The MARK statement must 


allocate enough space to store the entire program. 


Binary programs use more space on magnetic tape than ASCII programs. The amount 
depends on the size of the file. Because the SPACE function allocates the approximate 
space for an ASCII program, the SPACE function may or may not allocate enough space 
to hold the entire program in binary. If you try to execute the BSAVE routine and the 
selected file is not large enough to hold the current program, error message number 48 is 
displayed on the screen. To make sure enough space is allocated to hold the current pro- 


gram, you can use one of the following methods: 


METHOD 1 Allocate all of the space the program uses in RAM. For example: 
MARK 1, 32000 - MEM 
The MEM function returns the number of bytes still available 
in RAM. By subtracting this amount from the total storage 
capacity of RAM, the remainder is the amount of space the 
current program occupies in RAM. This remainder is also the 


amount of space needed to store the program on tape. 


METHOD 2 Allocate enough space to hold the entire contents of RAM. 
For example: 
MARK 1,32000 
This method may waste a lot of space on the tape if the current 


program is small. 
NOTE 
The actual storage capacity of RAM is approximately 32K or 


32000 bytes. Of this 32000 bytes, about 2000 are reserved 


by the microprocessor for a work area. The MEM function 
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BSAVE ROUTINE 


Syntax Form: 


[Line number] CALL ‘ BSAVE \ (,1/O address) 


string variable 


Descriptive Form: 


[Line number] CALL routine name [,1/O address] 





PURPOSE 


The BSAVE (Binary SAVE) routine stores the current BASIC program on the specified 
output device in binary code. If an output device is not specified, the internal magnetic 


tape unit is chosen by default. 


EXPLANATION 


The Random Access Memory (RAM) holds the current BASIC program in binary code. 
When executed, the BSAVE routine sends a copy of the current program to the output 
device in binary code. Like the SAVE statement, BSAVE does not alter assigned values 
of variables or system environmental conditions. The CALL “BSAVE” statement can be a 
step in the program being stored, or it can be executed directly from the Graphic System 


keyboard. 
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These statements store the current program as a secret program in binary code in file 1. 


From a TLIST command, you can see that the file header is marked BINARY PROGRAM 
SECRET as shown below. 


TLIST . 
1 BINARY PROGRAM SECRET 1624 
2 LAST 768 


To mark the header BINARY PROGRAM SECRET when using an external tape drive, 
such as the 4924, you must specify the |/O address of the tape drive ina SECRET com- 
mand. This allows the external device to recognize the program as secret and mark the 


header accordingly. 


For example, if the 4924 tape drive is assigned to device number 2 on the General Purpose 


Interface Bus, the following program stores the current program as secret in binary code. 


SECRET @ 2: 
SECRET 

FIND @ 2:1 

CALL “BSAVE"’,2 


USING AUTO LOAD WITH BINARY FILES 


The AUTO LOAD key finds the first file on the magnetic tape and then executes an OLD 
command automatically. The OLD command loads and executes the first program on 
file. Because the OLD command does not work with binary program files, the first pro- 


gram on tape must be an ASCII! program, if you use the AUTO LOAD key. 
By storing an ASCII program that finds and loads a binary program in file 1, you can 
use the AUTO LOAD key to automatically find and load a binary program. The following 


example shows how this can be done. 


Using a TLIST command, you can see the programs stored on the magnetic tape. 
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only considers the 30000 bytes that are free. By using 32000 
as a constant, you can make sure that you consider all of the 


possible storage capacity of RAM. 


Using the internal magnetic tape unit by default, the following example stores the current 


BASIC program in binary code. 


FIND 0 

MARK 1,1000 
FIND 1 

CALL “BSAVE” 


Since this is the first file on the tape, the read/write head is positioned at the beginning 

of the tape by the FIND 0 statement. The MARK statement allocates space on the tape 
to hold 1000 bytes. The read/write head is then positioned to the start of the allocated 
space by the FIND 1 statement. The CALL ‘‘BSAVE” statement then transfers a binary 


copy of the current program to file 1. 


If you now execute a TLIST statement, you can verify the action performed. Note the 
header of the internal tape for file 1 is marked BINARY PROGRAM. 


TLIST 
1 BINARY PROGRAM 1024 
2 LAST 768 


USING SECRET WITH BINARY PROGRAMS 


The current BASIC program can be marked SECRET and also be stored in binary code. 
The rules for using SECRET are the same as with any ASCII file. For example: 


SECRET 


FIND 1 
CALL “BSAVE” 
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BOLD ROUTINE 


Syntax Form: 


“BOLD” 
string variable 


} [,1/O address] 


[Line number] CALL 
Descriptive Form: 


[Line number] CALL routine name [, |/O address} 





PURPOSE 


The BOLD (Binary OLD) routine loads a copy of the program stored in binary code into 
the Random Access Memory (RAM) from the specified input source. |f an input source 


is not specified, the internal magnetic tape unit is chosen by default. 


EXPLANATION 


Any program that is in the machine dependent binary code of the 4051 Graphic System 
can be loaded into RAM by using the BOLD routine. A program stored by the BSAVE 
routine is usually the program retrieved by the BOLD routine. BOLD also loads BINARY 


SECRET programs, and the retrieved programs remain secret. 
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Cal PD 


HSCIT PROGRAM 
BINARY FROGRAM 
LAST 


File 1, an ASCII program file, contains this program: 


100 FIND 2 
110 CALL “BOLD” 


When the AUTO LOAD key is pressed, the program in file 1 is loaded and executed. 


This program in turn finds and loads the desired binary program in file 2. Since the 


CALL “BOLD” statement is executed under program control, the binary program in file 


2 is executed after being loaded. 


SPECIFYING AN OUTPUT DEVICE 


The current program can be stored in binary code on any peripheral device in the system 


by specifying the desired primary address in the CALL “BSAVE” statement. For 


example: 


290 FIND @ 19:5 
300 CALL “BSAVE"’’,19 


This statement sends the current program in binary code to file number 5 on device 19 


on the General Purpose Interface Bus. 


i 





et, 


CAUTION 


Fa a a 





Unless the internal tape file is closed when executing a BSAVE to 


an external device, data may be inadvertently written on the 


internal magnetic tape. Ensure the internal tape file is closed 


before executing a BSAVE to an external device. 


REV. A, JAN. 1977 
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BAPPEN ROUTINE 


Syntax Form: 


“BAPPEN”, 


string variable, (1/0 address;] tine number [,increment] 


[Line number] CALL 


Descriptive Form: 


target line number line number | 


[Line number] CALL routine name, [1/O address;] . 
in current program —,increment 





PURPOSE 


The BAPPEN (Binary APPEND) routine inputs BASIC statements stored in binary code on 
the specified peripheral device and attaches those statements to the program currently in 


RAM. The internal magnetic tape unit is selected by default if an input device is not 


specified. 


EXPLANATION 


The BAPPEN routine follows the same procedure as the APPEND statement. First, the 
read/write head of the input device must be positioned to the beginning of a binary pro- 
gram file by using the FIND command. Then, when the BAPPEN routine is executed, the 
given target statement is overwritten by the first statement coming from the peripheral 
device. The newly appended statements and any statements that originally followed the 
target statement are renumbered, stating with the target line number. The line numbers 


are incremented by either the given increment, or by the default of 10 if an increment ts 


not specified. 
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LOADING A BINARY PROGRAM FROM MAGNETIC TAPE 


To retrieve a binary program from magnetic tape, the read/write head of the output device 
must first be positioned at the beginning of a binary program file. For example, using the 


internal magnetic tape unit by default, the statements: 


FIND 1 
CALL “BOLD” 


load the binary program from file 1 into RAM. The FIND statement locates the beginning 
of file 1. The BOLD routine erases everything currently in RAM, then transfers the binary 
file into RAM. The loaded program is ready to be executed by a RUN statement or edited. 
Like the OLD command, if the BOLD routine is executed under program control, a RUN 


statement is automatically executed after the program is loaded into RAM. 


SPECIFYING AN INPUT DEVICE 
Any peripheral device holding a binary program can be specified as the input source for the 
BOLD routine by designating the appropriate 1/O address. The following statement 


specifies file number 6 on device number 2 on the General Purpose Interface Bus. 


290 FIND @ 2:6 
300 CALL "“BOLD”,2 
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LINK ROUTINE 


Syntax Form: 


“EINK”, 


[Line number] CALL es variable, 


\ {1/O address:] line number 


Descriptive Form: 


[Line number] CALL routine name, [I/O address;] line number of entry point 





PURPOSE 


The LINK routine transfers a stored binary program to RAM from the specified peripheral 
device without disturbing variables and associated values stored in RAM. If an input device 


is not selected, the internal magnetic tape unit is selected by default. 


EXPLANATION 


The LINK routine erases the program currently in RAM, but retains all variables and their 
currently assigned values. Then, a binary program is loaded into RAM from the specified 
input device. The BASIC line counter is set to the starting line number and execution 


starts at that point. 


These features allow you to break a large BASIC program into subroutines, store them on 
tape, and then structure the flow of execution by using the LINK routine. By using the 

LINK routine, the size of a program is no longer limited by the storage capacity of RAM. 
\f RAM cannot hold an entire program, you can break the program into sections and use 


the LINK routine to load and execute the program sections. 
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For example: 


200 REM APPEND STATEMENTS FROM FILE 2 HERE 
FIND 2 
CALL “*BAPPEN’’,200 


When these statements are executed, the FIND statement positions the read/write head 

at the beginning of File 2. File 2 must be a binary program file. (Error message number 
55 is printed if file 2 is not binary). The BAPPEN routine replaces line 200 with the first 
statement in file 2. The remaining statements in file 2 are then appended with line num- 


bers incremented by 10. 


BAPPEN WITH SECRET PROGRAMS 
lf the program you want to append is stored as a SECRET BINARY program on the periph- 


eral device, after executing the BAPPEN routine, the entire contents of RAM becomes 
secret. This includes all original statements and the statements appended. 


SPECIFYING AN INPUT DEVICE 


Any peripheral device storing a binary program can be specified as the input device. The 


device is specified like this: 


FIND @ 15:3 
CALL ‘‘BAPPEN”, 15;650,50 


This statement selects file number 3 on device number 15 on the General Purpose Inter- 


face Bus as the input device. Line numer 650 is the target statement in the current pro- 


gram stored in RAM. The renumber increment is 50. 
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FIND 1 

CALL “BOLD” 

LIST 

100 INIT 

110 AS="CHANGE” 

120 PRINT “SOME THINGS NEVER ”;A$ 


Now, the program is executed. 


RUN 
SOME THINGS NEVER CHANGE 


This binary program is used in the following examples by the LINK routine. 


100 A$="STAY THE SAME” 
110 FIND 1 
120 CALL “LINK”, 120 


In line 100, A$ is given the value “STAY THE SAME”. The binary file is then located by 
the FIND statement. The LINK routine directs execution to begin in fine 120 of the binary 


program. So... 


RUN 
SOME THINGS NEVER STAY THE SAME 


The stored binary program is now in memory. No statements have been changed. A$ 
retains the value “STAY THE SAME” because execution started in line 120. 


Notice the INIT statement in line 100 of the binary program. If the LINK routine is 


specified by 


CALL “LINK’’,100 
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LOCATING THE LINK POINT 


To use the LINK routine, the read/write head of the peripheral device must first be posi- 


tioned at the beginning of a binary program file. 


To locate the beginning of a binary program file, use the FIND statement. After the LINK 
routine is executed, variables and values are kept and the current program is erased. The 


entire binary program (previously located by the FIND statement) is then loaded into RAM. 
if the LINK routine is executed in immediate mode, execution stops after the program is 


transferred. The loaded program can then be executed by entering a RUN statement and 


pressing the RETURN key. Execution begins with the starting line number. 


if the LINK routine is executed under program control, then after the binary program is 


transferred, execution automatically begins with the starting line number. 


USING LINK WITH SECRET BINARY PROGRAMS 


If a LINK routine brings a SECRET BINARY program into RAM, then the entire contents 


of RAM are made secret. As always, secret programs can only be executed. 


SAMPLE PROGRAMS 


To use the LINK routine a BASIC program must be stored in binary code. As an example, 
the following program is stored in file 1 on the internal magnetic tape unit. Since this file 
is stored in binary code, the file is retrieved by using the BOLD routine. A LIST statement 


is executed to show the program. 
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SUMMARY 


1. The BINARY PROGRAM LOADER adds four routines to the Tektronix 4051 Graphic 
System. Three of the routines (BSAVE, BOLD, BAPPEN) perform the same basic 
Operation on binary files as the commands (SAVE, OLD, APPEND) perform on 
ASCII! program files. The following chart sumarizes the four routines and their 


relationships to other statements. 





BINARY ASCII 
BSAVE SAVE 
BOLD OLD 
BAPPEN APPEND 
LINK 


2. All routines are specified by using a CALL statement. 

3. All routines access binary program files. 

4. The read/write head of the peripheral device must be positioned to the beginning 
of the binary program file before executing a BINARY PROGRAM LOADER 


routine. Use the FIND statement to locate the beginning of a file. 


5. Binary programs may be also declared SECRET. The SECRET BINARY program 
remains secret when used by a BINARY PROGRAM LOADER routine. 
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then execution would begin with line 100. The INIT statement puts al! variables in an 
undefined state, including any variables retained by the LINK routine. A$ would be 


undefined. Then in line 110 of the binary program, A$ would be set to “CHANGE”. 


If the LINK routine ts specified by 


CALL “LINK’’,110 


then A$ is retained with the value “STAY THE SAME”, but is set to “CHANGE” in line 
110 of the binary program. 


From these exampies, you can see how execution and program flow can be controtled by 


using the LINK routine. 


SPECIFYING THE INPUT DEVICE 


The LINK routine uses any peripheral device in the Graphic System. You can specify 
the desired peripheral device like this: 


FIND @ 2:8 
CALL “LINK”, 2;150 


In this case, file number 8 on device number 2 is selected as the input device on the Gen- 
eral Purpose Interface Bus. Execution begins with line 150. 
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MESSAGE 
NUMBER 


14 


15 


32 


38 


Appendix B 


ERROR MESSAGES 


ERROR MESSAGE 
———— PROBABLE CAUSE AND SOLUTION 


INVALID COMMAND ARGUMENT 

———— A parameter in the CALL statement for the BAPPEN 
routine is invalid. Check the syntax of the CALL state- 
ment for the BAPPEN routine. 


INVALID COMMAND ARGUMENT 
———— The line number specified in the CALL statement for the 
BAPPEN routine does not exist. List the program to find 


the correct line number. 


CALL NAME INVALID 

——-—-— The routine name is misspelled. Check the spelling of 
the selected routine name. 

———— Asstring variable appears in a CALL statement with a 
value of no recognizable routine name. Check the spelling 
of the value assigned to the string variable. 

———— The 4051R05 ROM pack is not installed in the machine. 


See the installation instructions included in this manual. 


PROGRAM IS SECRET 
———— You attempted to modify or output a secret program. 


Secret programs can only be executed after being retrieved. 
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61 INVALID OPERATION ON AN OPEN FILE 
———— A MARK statement was executed and the read/write 
head of the peripheral device was not positioned at the 
beginning of a file. Use the FIND statement to position 


the read/write head. 
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48 


51 


52 


55 


57 


EOF ON UNIT O IN IMMEDIATE LINE 
———— You tried to store a program in a file that is not large 
enough to hold the entire program. Either allocate more 


space for that file or locate or make a file with more space. 


INVALID LINE NUMBER 
———— The line number specified cannot be found or is invalid. 


List the program to find the correct line number. 


MAG TAPE FILE NOT FOUND 
———— Either the specified magnetic tape file does not exist or 
an attempt has just been make to KILL the LAST (dummy) 
file. 
Issue a TLIST to find the desired file number. 


MAG TAPE ERROR 

———— The read/write head is not at the beginning of a file before 
accessing that file. Execute a FIND statement to position 
the read/write head before using the file. 

-—~-—— You tried to use an ASCII program file with one of the 
BINARY PROGRAM LOADER routines. The BINARY 
PROGRAM LOADER routines only operate on a BINARY 
PROGRAM file. 

——-—-— You tried to use an ASCII program file command with a 
BINARY PROGRAM file. BINARY PROGRAM files can 


only be used by the BINARY PROGRAM LOADER routines. 


MAG TAPE CARTRIDGE REQUIRED 
——-—-— A magnetic tape cartridge is not inserted into the magnetic 
tape unit. 


Insert a magnetic tape into the tape unit. 
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